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Abstract 

An explicit formula for the mean spectral measure of a random Jacobi 
matrix is derived. The matrix may be regarded as the limit of Gaussian 
beta ensemble (G/3E) matrices as the matrix size N tends to infinity with 
the constraint that Np is a constant. 
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1 Introduction 

The paper studies spectral measures of 
the form 

/A/’(0, 1) X2a 
J„= X2a AA(0,1) 

where the diagonal is an i.i.d. (independent identically distributed) sequence of 
standard Gaussian A/'(0,1) random variables, the off diagonal is also an i.i.d. se¬ 
quence of X 2 a-distributed random variables. Here X 2 a = X 2 alV^ with X 2 a 
denoting the chi distribution with 2 q; degree of freedom. As explained later, 
is regarded as the limit of Gaussian beta ensembles (G/3E for short) as the ma¬ 
trix size N tends to infinity and the parameter j3 also varies with the constraint 
that Nj3 = 2a. 

Let us explain some terminologies and introduce main results of the paper. 
A (semi-infinite) Jacobi matrix is a symmetric tridiagonal matrix of the form 

( ai bi 
hi 02 


^ , where o^ G M, 5^ > 0. 


'andom (symmetric) Jacobi matrices of 


X2a , {a > 0), 


1 



For a Jacobi matrix J, there is a probability measure /r on K such that 

[ x'^dn = (ei) = J^{1, 1), fc = 0 , 1, ..., 

Jr 

where ei = G Here {u,v) denotes the inner product of u and 

V in while (p,, f) ■= J fd^ will be used to denote the integral of a function 
/ with respect to a measure /r. Then the measure fj, is unique if and only 
if J, as a symmetric operator defined on Dq = {x = (xi, 3:2, ■ ■ • : Xk = 

0 for k sufficiently large}, is essentially self-adjoint, that is, J has a unique self- 
adjoint extension in When the measure /i is unique, it is called the spectral 
measure of J, or more precisely, the spectral measure of (J, Ci). It is known 
that the condition 


1 


Er = 


00 


i=l 


implies the essential self-adjointness of J, [3 Corollary 3.8.9]. 

For the random Jacobi matrix J^, the above condition holds almost surely 
because its off diagonal elements are positive i.i.d. random variables. Thus 
spectral measures /ia are uniquely determined by the following relations 

(Ma,a;'=) = J^(l,l), fc = 0,l,.... 

Then the mean spectral measure p,a is defined to be a probability measure 
satisfying 

ip-aj) =E[(/ia,/)], 

for all bounded continuous functions / on K. It then follows that 

(fla.X^) =E[{fXa,X^)], A: = 0,1,..., 

provided that the right hand side of the above equation is finite for all k. 

The purpose of this paper is to identify the mean spectral measure fia ■ Our 
main results are as follows. 

Theorem 1. (i) The mean spectral measure pLa coincides with the spectral 

measure of the non-random Jacobi matrix Aa, where 

0 Va J-1 

^ _ I \/o J- 1 0 \/q! J- 2 


(ii) The measure pa has the following density function 

e-^V2 1 


= 


^ l/a(y)P’ 


where 


Uiy) = 


fait)e"y*dt, fa{t) = T^Jwr^t' 


r(Q!) 
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Let us sketch out main ideas for the proof of the above theorem. To show 
the first statement, the key idea is to regard the Jacobi matrix Jq, as the limit 
of G/3E as the matrix size N tends to infinity with iV/3 = 2a. More specifically, 
let Tjq{j3) be a finite random Jacobi matrix whose components are (up to the 
symmetry constraints) independent and are distributed as 


Tn{P) 


/7V'(0,1) X(N-i)p 
X(W-l)/3 AA(0, 1) X{N-2)P 


\ 


V 


X/3 AA(0,1)/ 


Then it is well known in random matrix theory that the eigenvalues of Tjv(/3) 
are distributed as G/3E, namely. 


N 

(Ai,..., Aat) oc \Xk — Xj\^. 

1=1 l<j<k<N 


Moreover, by letting N —>■ oo with /3 = 2a/N, the matrices Tn{13) converge, in 
some sense, to Jq. That crucial observation together with a result on moments 
of G/3E ([^ Theorem 2.8]) makes it possible to show that p,a coincides with the 
spectral measure of Aa- 

The next step is to establish the following self-convolutive recurrence for 
even moments of p,a, 

n—1 

Un{c^) = (2n - l)un-i{a) + a ^ Ui{a)un-i-i{a), 

where Un{a) is the 2nth moment of /Jq. Note that its odd moments are all 
vanishing because the spectral measure of Aa is symmetric. Finally, the explicit 
formula for Jla is derived by using the method in [4j. 

The paper is organized as follows. In the next section, we mention some 
known results on G/3E needed in this paper. In Section 3, we introduce the 
matrix model and step by step, prove the main theorem. 


2 A result on Gaussian /3-ensembles 

The Jacobi matrix model for G/3E, a finite random Jacobi matrix, was dis¬ 
covered by Dumitriu and Edelman First of all, let us mention some pre¬ 
liminary facts about finite Jacobi matrices. Assume that J is a finite Jacobi 
matrix of order N (with the requirement that the off diagonal elements are 
positive). Then the matrix J has exactly N distinct eigenvalues Ai, A 2 ,..., Ajy. 
Let vi,V 2 , ■ ■ ■ ,vn be the corresponding eigenvectors which are chosen to be an 
orthonormal basis in . Then the spectral measure /r, which is well defined 
by (/i, a;^) = J^{1, 1), fc = 0,1,..., can be expressed as 

N 

M = qj = \vj{l)\, 

i=i 
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where 5 a denotes the Dirac measure. It is known that a finite Jacobi matrix of 
order N is one-to-one correspondence with a probability measure supported on 
N points, or a set of Jacobi matrix parameters {ai}fLi, {bj}^Si is one-to-one 
correspondence with the spectral data {Xi}iLiA<lj}^=i- 

The Jacobi matrix model for G/?E is defined as follows. Let {ai}fLi be an 
i.i.d. sequence of standard Gaussian A/’(0,1) random variables and be 

a sequence of independent random variables having x distributions with param¬ 
eters {N — 1)13, {N — 2)13 ,..., 1, respectively, which is independent of 
Here Xk, for k > 0, denotes the distribution with the following probability den¬ 
sity function 


2 

rW^ 


k — 1 —u‘^ 

u e 


,u> 0, 


which is nothing but Xfc/v^, or the square root of the gamma distribution with 
parameter (fc/2,1). We form a random Jacobi matrix Tm{13) from {ai}fLi and 
as follows, 


TnW) 


/A/’(0, 1) X(Af-l)/3 

X(Ar-l)/3 -^5(0,1) X{N-2)I3 


\ 


V 


X /3 A5(0, 1)J 


Then the eigenvalues and the weights {qj}^^i are independent, with 

the distribution of the former given by 

N 

(Ai, A 2 ,..., Aat) oc |Afe —Ajl^, 

l<j<k<N 


and the distribution of the latter given by 


{qi,q2,...,qN) oc 


qN 


N 

n 




N 

{q, > = 1 ). 

i=l 


It is also known that q = {qi,..., q^) is distributed as a vector (x/ 3 , ■. ■ ,Xi 3 ) 
with i.i.d. components, normalized to unit length. 

The trace of Tjq{f3)'^ and r/v(/3)"(l, 1) can be expressed in term of the spec¬ 
tral data as 


N N 

Tr(T^(/3)") = ^ A,", 1) = ^ 9 ^,". 


Consequently, 


N 


N 


E[T^(/3)"(l, 1)] = E[^ g|A,"] = ^ E[< 7 |]E[A,"] = ^ ^ E[A,"] 
1=1 i=i i=i 

= lE[Tr(X^(/3)”)]. 


N 
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In the rest of this section, for convenience, we use the parameter /3 = /3/2. 
Let mp{N,P) = E[r/v(2/3)^P(l, 1)]. It is clear that mp{N,P) is a polynomial of 
degree p in N, and thus mp{N,f3) is defined for all G K. Then a result for 
the trace of Tjv(/3)"' can be rewritten for mp{N, fi) as follows. 

Theorem 2 (cf. [5J Theorem 2.8] and [71 Theorem 2]). It holds that 
mp{Nj) = i-irpPmpi-mJ-^). 


Observe that j3 ^mp{N,j3) is the expectation of the 2pth moment of the 
spectral measure of the following Jacobi matrix 


^r^(2/3) = - 

V/3 


As ,5 —>■ oo, it holds that 


^ A/’(0, 1) X(Ar-i)2/3 ^ 

X(Ar_i)2/3 1) X(N- 2 ) 2 $ 


V 


aa( 0 ,i)_q Xk20 _ / r(fc^,i) 


1/2 


X20 Af(0,1)/ 


Vk ( in L'^ for any g > 1). 


The convergences also hold almost surely. Therefore as /3 ^ oo, 


— j=TNi‘2P) —> 


/ 0 


0 




\ 

=: Hjm. 


1 0 / 


Here the convergence of matrices means the convergence (in L'^) of their ele¬ 
ments. Let hp{N) = iL^(l,l) for N > p. Then hp{N) is a polynomial of 
degree p in TV so that hp{N) is defined for all N G M.. The above convergence 
of matrices implies that for fixed p and fixed N, 


Let 


hp{N) = lim $-Pmp{Nj). 
^—>•00 


-^a 


/ 0 \/Fn: 

v/FFT 0 



( 1 ) 


and let Up{a) = 1). Then Up{a) is also a polynomial of degree p in a. In 

addition, it is easy to see that 

Up{a) = {-l)Php{-a). (2) 

As a direct consequence of Theorem and relations Q and ([^, we get the 
following result. 
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Proposition 3. As N ^ oo with /3 = /3(-/V) = a/N, 

mp{Nj) Up{a) = 1 ). 


3 Random Jacobi matrices related to Gaussian 
13 ensembles 

3.1 A matrix model and proof of Theorem [^i) 

Consider the following random Jacobi matrix 

/7V(0,1) X2a \ 

J„= X2a Af(0,l) X2a 


where all components are independent random variables. More precisely, the 
diagonal {ai}“ i is an i.i.d. sequence of standard Gaussian A/’(0,1) random vari¬ 
ables and the off diagonal {bj}^i is another i.i.d. sequence of X 2 a random vari¬ 
ables. Then the spectral measure fXa of Jq, exists and is unique almost surely 
because 


E 




oo(almost surely). 


The mean spectral measure Ha is defined to be a probability measure satis¬ 
fying 




for all bounded continuous functions / on M. Then Theorem [^i) states that 
the measure Jla coincides with the spectral measure of {Aa, ei). 


Proof of Theorem [^i) . Note that the spectral measure of A ^, 
sure /r satisfying 




a probability mea¬ 


ts unique because 


Also, it is clear that 


E 


i=i 


1 

VoTj 


= CX). 


(/i„,x'') = E[(/rQ,,a:'')], A: = 0,1,..., 


because E[(/ra, |a;|^)] < oo for all fc = 0,1,.... Therefore, our task is now to 
show that for all fc = 0,1 ,..., 


(a1„,x'=)=A^(1,1). 


(3) 
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We consider the case of even k first. For any fixed j, all moments of the 
X(N-j)2$ distribution converge to those of the X2a distribution as iV —?> oo with 

f3 = a/N. Thus for fixed p, as N —>■ oo with /3 = a/N, 

mp(iV,/3) = E[Tjv(2/3)2p(l, 1)] ^ E[J2P(1,1)] = E[(/r„,a;2P)]. 
Consequently, for even k, namely, k = 2p, 

(Aa,x'=)=A^(l,l), 
by taking into account Proposition 

For odd k, both sides of the equation ([^ are zeros. Indeed, = 0 

when k is odd because the diagonal of is zero. Also all odd moments of Pa 
are vanishing, 

= E[(/ia, = 0 , 

because the expectation of odd moments of any diagonal element of are zero. 
The proof is completed. □ 


3.2 Moments of the spectral measure of 

Recall that 

M„(a) = A2"(l,l),n = 0 ,l,.... 

Proposition 4. (i) Un{a) is a polynomial of degree n in a and satisfies the 

following relations 


Un{a) = (a + 1) u^{a + l)un-i-z{a), n > 1, 

uoia) = 1 . 


(4) 


(ii) {'«n(<a)}^o satisfies the following relations 

\un{a) = {f2n-l)un-i{a)v,i{a)un-i-i{a), n > 1, 

[uo(a) = 1 . 

Remark 5. The sequences {««(«)}n> 0 j for a = 1 and a = 2, are the sequences 
A000698 and A167872 in the On-line Encyclopedia of Integer Sequences [5], 
respectively. Relations @ and © as well as many interesting properties for 
those sequences can be found in the above reference. In the proof below, we 
give another explanation of u„(a) as the total sum of weighted Dyck paths of 
length 2n. 

Proof. In this proof, for convenience, let the index of the matrix Aa start from 
0. Since the diagonal of Aa is zero, it follows that 

2n-l 

^a”(0,0)= ^ l[Aa{t„t, + A, 

{io,ii,...,i2Ti}eS2n j=0 
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where S 2 n denotes the set of indices {zq, ii,..., i 2 n} satisfying that 

Zq 0 , Z 2 n 0 , ij ^ 0 , 

l*i+i - = 1>J = 0,1,.. .,2rz- 1. 

Each element in S) 2 n corresponds to a path of length 2n consisting of rise steps 
or rises and fall steps or falls which starts at ( 0 , 0 ) and ends at ( 2 rz, 0 ), and stays 
above the x-axis, called a Dyck path. We also use D 2 n to denote the set of all 
Dyck paths of length 2n. 

A Dyck path p is assigned a weight w{p) as follows. We assign a weight 
(a + fc + 1) for each rise step from level k to k + 1, and the weight w{p) is the 
product of all those weights. Then 

Unia) = Al^{0,0) = w{p). 

PeS2n 



Figure 1: A Dyck path p with weight w{p) = (a + 1 )^(q; + 2)^{a + 3)(q; + 4). 

Let SJn be the set of all Dyck paths of length 2n which do not meet the 
x-axis except the starting and the ending points. Let 

Vn{a) = Y 

Since each Dyck path p = (zq, zi,..., Z 2 n-i, * 2 n) € ® 2 n is one-to-one correspon¬ 
dence with a Dyck path g = (zi — 1, Z 2 — 1,..., Z 2 n-i — 1) of length 2(n — 1), it 
follows that 

Vnia) = (a -I- l)zt„_i(a -I- 1). 

Moreover, let 2z be the first time that the Dyck path p meets the x-axis. Then 
either z = zz or the Dyck path p is the concatenation of a Dyck path in (1 < 






i < n), and another Dyck path of length 2(n — i). Thus, 


Un{a) = Vn{a) + ^ Vi{a)un-iia) 


= (a + 1 )'U„_i(q; + 1) + ^(a + + l)u„_i(a) 

n—1 

= (a + 1) ^ u^(a + l)'Un-i-i(a). 

i^O 

The proof of (i) is complete. We will prove the second statement after the next 
lemma. □ 

Lemma 6. Let a > 0 be fixed. Let {a„} he a sequence defined recursively by 

f — (2?T, T Oi Q CliCln—1 — i: n ^ 1, 

\ao = 1 . 

Let {bn} be a sequence defined by the following relations &o = 1, 

n—1 

«n = (a + 1) ^ bian-i-i, n > 1. (7) 

i^O 

Then {bn} satisfies an analogous recursive relation as {a„}, 

ibn = i2n-l)bn-i + + bibn-i-i, n > 1, 

\6o = l. 

Proof. Consider the field of formal Laurent series over M, denoted by IR((X)), 
M((X)) = < f{X) = ^ c„X” : c„ G K, c„ = 0 forn < no > . 

I nGZ J 

The addition is defined as usual and the multiplication is well defined as 

nGZ \ieZ / 

for f{X) = J2cnX'^,9iX) = J2dnX^ G M((X)). The quotient f{X)/g{X) is 
understood as f{X)g{X)~^ for g{X) 0. The formal derivative is also defined 
as 

/'(X) = ^c„nX"-iGM((X)). 


Now let 
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It is straightforward to show that the recursive relation ([^ is equivalent to the 
following equation 

f{X) - 1 = 2X^f{X) + Xf(X) + aXfiX). 

In addition, the relation 0 leads to 

fix)-i 

^ {a + l)Xf{X)- 

Finally, we can easily check that g{X) satisfies 

giX) - 1 = 2X^g'iX) + Xg{X) + (a + l)Xg\X), 

which is equivalent to the recursive relation Q. The proof is complete. □ 

Proof of Proposition^ii). When a = 0, it is well known that m„( 0) is the 2nth 
moment of the standard Gaussian distribution, and is given by 

M„(0) = ( 2 n-I)!!. 

Consequently, the conditions in Lemma are satisfied for a„ = Un(0),bn = 
Un(l) and a = 0. It follows that the recursive relation then holds for 
a = I. Continue this way, it follows that the recursive relation 0 holds for any 
a G N. We conclude that it holds for all a because of the fact that {«„(«)} is 
a polynomial of degree n in a. The proof is complete. □ 


3.3 Explicit formula for the spectral measure of Aa, proof 
of Theorem [^ii) 

In this section, by using the method of Martin and Kearney [1], we derive the 
explicit formula for the mean spectral measure fia from the relation ([^, 

j Un{a) = (2n - 1 )u„_i(q!) + Ui{a)un-i-i{a), n > 1, 

[iio(a) = I. 


Recall that it„(a) = and fia is a symmetric probability measure. 

Let us extract here the main result of [4]. The problem is to find a function 
V for which 

x^~^i'{x)dx = Un, n = 1 , 2 ,..., 
where the sequence {u„} is given by a general self-convolutive recurrence 
Un — (oiTZ T 02 )^ 71 —1 T 0^3 ^ ^ 2, 

Ui = I, 



oi, 02 and 03 being constants. Then the solution is given by (Eq. (13)-Eq. 

in ID), 

k{kx)~^e~^^ I 

r(a + l)f (a — 6+1) Uii{kxY + Ui{kx)'^ ’ 


(16) 
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where, 


- (cos7r6) ^^_^ - a; 2 - b-,x)\ 

r(a) J 

Ui{x) = (sin7r6)e~^ ^^^ — a; 2 — 6; a;), 

i (a) 


and k = 1/ai, a = a^/ai, b = —1 — a 2 /ai, provided ai ^ 0. Here ii^i(a; b; z) is 
the Kummer function. 

The sequence {un(a)}n>o is a particular case of the self-convolutive recur¬ 
rence with parameters a\ = 2, a 2 = ~3 and = ol. Note that our sequence 
{M„(a)} starts from n = 0, and thus a 2 = —3. By direct calculation, we 
get k = 1/2, a = a/2, and b = 1/2. Therefore, the function Vaix) for which 
Un{<x) = /q°° x'^dva{x)dx, u = 0, 1 ,..., is given by 


Va{x) 

where 


^/2V{% + l)r(f + \)y/^ Ur{x/ 2Y + Ui{x/2.f 


r(l) 


1 a 1 


Ur{x) e 2’2’*^’ 

Ui{x) = ^;a:). 


r(f) 


2 2 


a: > 0, 


( 10 ) 

( 11 ) 


It is clear that v^ix) > 0 for any a; > 0. Now it is easy to check that the function 
ftaiy) defined by 

P-aiy) = lyWaiy"^), y g K, 

satishes the following relations 


y y-aiy)dy = 0, 



fJ-a{y)dy 


Un{a), n = 0,1,..., 


In other words, p-aiy) is the density of the mean spectral measure /Iq, with 
respect to the Lebesgue measure. 

We are now in a position to simplify the explicit formula of fla- Let 


VR{y) = 


r(f+ i)r(f+ i)V^' 


nh) 


= 2 2r(a-|-l)2 

Vi{y) = - 


V 

r(i) 


Vr{,V^/2), 


r(f+ i) 
r(f+ i)r(f+ i)V^' 


_f ^ ,1 a 1 . 

XyS ^ iTi(-— —; 


r(i) 


ti/(yV2) 


= -2 ? 3r(a + l)5^J^ye iFi(I-|;5;|-). 


II 













Here, in the above expressions, we have used the following relation for Gamma 
function 


r(i) 


Then fJ-aiy) can be written as 


y-aiy) = 


e 2 


Vniy)^ + Viiyr' 


Next, we will show that Vii{y) and Vi{y) are the Fourier cosine transform 
and Fourier sine transform of 


fa{t) = TT 


-t' 


e 2 


r(a) v^’ 


respectively. Let us now give definitions of Fourier transforms. The Fourier 
transform of a function /: K —>■ C is defined to be 


1 r°° 

J^{f){y) = / f{t)e^y*dt, {y G 

V j —OO 


and the Fourier cosine transform, the Fourier sine transform are defined to be 

/2 roo 


^cU){y) = 


f{t)cos{yt)dt, (y > 0), 


^s{f){y) = ^J^ f{t)sm{yt)dt, {y>0), 

respectively. Then those transforms are related as follows 

f ^if)iy) = ^c{f){y), {y > 0), if f{t) is even, 
\j"if)iy) = 'i-J"s{f){y), iy>0), if/(t) is odd. 

For a > 0, we have (cf. Formula 3.952(8) in [3]) 

Then by some simple calculations, we arrive at the following relation 
Vniy) = ^c{fa{t)){y), y>0. 

Similarly, 

Viiy) = Tsifa{t)){y), y>0, 
by using Formula 3.952(7) in [3], 
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By definitions, Vii{y) is an even function and Vi{y) is an odd function. Thus 
the following expression holds for all y gM., 

VR{y) + iVi{y) = ^ J fa{t){cos{yt) +ism{yt)dt 

= \li Ut)e^^dt=:Uy)- 

Consequently, 

VR{yf + Vj{yf = \Uv)?, 

which completes the proof of Theorem Bii)- 

We plot the graph of the density {y) for several values a as in the following 
figure by using Mathematica. It follows from the Jacobi matrix form that the 
spectral measure of converges weakly to the semicircle law as a tends 

to infinity. Note that the semicircle law, the probability measure supported on 
[—2, 2] with the density 

^\/4-a:2, (-2 < a; < 2), 

is the spectral measure of the following Jacobi matrix 

/O 1 \ 

1 0 1 
10 1 

V ■■•/ 



Figure 2: The density y-a{y) for several values a. 
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Remark 7. When a in a positive integer number, we can give even more explicit 
expressions for VR{y) and Vi{y). 

(i) a = 2n,n € N. In this case, fa{t) is an odd function. Therefore 

Vi{y)=Ts{Ut)) = -iJ^{Ut)). 


Note that 


J-(e-T) = e- 


Therefore, for integer a > 1, 


Ja—1 




Consequently, 


Vi{y) = -i°‘Tr 


a (i“ ^ , _id 1 
-(e 2 ) 


r(a) dy°^ 


-1 




a !d(i“^,_yd.e 2 

-e 2 ^--(e 2 )- 

y — 1 ^ ' 


r(a) dy° 




a 6 2 

-Hea-i- 


Here Hem denotes probabilists’ Hermite polynomials. 

(ii) a = 2n + 1. This case is very similar. Since fa(t) is an even function, it 
follows that 


Vniy) = MUt))iy) = HUt)) = 


r(a) 


H C(x— 1 


e 2 

V^' 
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